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Abstract 

Following recent results from the SNO solar neutrino experiment and the K2K 
long-baseline neutrino experiment, the combined existing data on neutrino oscil- 
lations now point strongly to a specific form for the lepton mixing matrix, with 
effective bimaximal mixing of and v T at the atmospheric scale and effective tri- 
maximal mixing for v e with f„ and v T at the solar scale (hence 'tri-bimaximal' 
mixing). We give simple mass-matrices leading to tri-bimaximal mixing, and 
discuss its relation to the Fritzsch-Xing democratic ansatz. 
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1. Motivation 



The first data from the Sudbury solar neutrino experiment (SNO) [|TJ have dramati- 
cally confirmed the long-standing HOMESTAKE solar neutrino result with regard 
to the high-energy charged-current z/ e -rate (SNO/BP2001 = 0.347 ± 0.029, HOME- 
STAKE/BP2001= 0.337 ± 0.030). At the same time, the comparison with the exist- 
ing rather precise SUPER-K (SK) result for solar-neutrino electron elastic scattering 
(SK/BP2001 = 0.459 ±0.007 0]) (which includes a neutral-current contribution) has 
provided a significant first cross-check of the Bahcall-Pinsonneault (BP2001) [|j] stan- 
dard solar model calculation of the 8 B flux in these experiments, so that it now seems 
reasonable to conclude that the suppression for 8 B neutrinos is probably close to 
~ 1/3. In detail, the SNO, HOMESTAKE and SUPER-K experiments have different 
thresholds and response functions (eg. HOMESTAKE is expected to include a ~ 15% 
contribution from 7 Be neutrinos) but such effects are readily taken into account 0]. 

By comparison, the low energy charged-current v e rate as sampled in the gallium- 
based experiments SAGE , GALLEX |7| and GNO [§] is known to be less suppressed 
(SAGE/BP2001 = 0.59 ±0.07, (GALLEX GNO)/BP2001 = 0.58 ±0.07). The gallium 
experiments are sensitive to neutrinos from the pp chain and are only marginally 
affected by 8 B neutrinos (with expected signal contributions of ~ 60% pp, ~ 30% 7 Be 
and ~ 10% 8 B in the standard model 0). We have previously emphasised, within the 
context of the original trimaximal mixing scenario |J, the consistency of the gallium 
suppression with 5/9 ~ 0.56 (this consistency survives today at the 1.2 a level, even 
allowing for the reduced 8 B contribution in gallium in the LMA solution, see below). 

Thus energy-dependence of the solar suppression is implicit, with the latest general 
fits || to the solar data favouring the so-called large-angle (LMA) MSW |10] solution. 
The long-standing small-angle (SMA) MSW solution is now essentially ruled out, while 
the so-called LOW and VO solutions are of marginal significance only ||. The LMA 
solution is illustrated in Fig. 1 for several possible mixing angles. In the LMA solution 
the base of the MSW 'bathtub' is arranged to account for the strong suppression at 
high-energy where matter effects dominate. At lower energies (for the same solar core 
density) the suppression reverts to its vacuum level outside the bathtub, accounting 
for the gallium data (the far high-energy end of the bathtub plays no role in the 
LMA solution). No significant day-night asymmetry is observed [|[] so that the latest 
LMA fits now prefer a mass-squared difference towards the higher end of the range 
~ 10 -5 — 10 -4 eV 2 (the curves of Fig. 1 are for a representative Am' 2 = 5 x 10~ 5 eV 2 ). 

Interestingly, the trimaximal model || is known [1 L] to predict a '5/9 — 1/3 — 5/9' 
LMA bathtub (see Fig. 3 of Ref. |L1|]) which could certainly be exploited to fit the 
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current solar data in isolation. However, within the trimaximal model, the associated 
mass-squared difference would then necessarily be the larger mass-squared difference 
(compare Figs. 2 and 3 of Ref. (JTI]]) and would thus be inconsistent with the data 
on atmospheric neutrinos, which seem to require a mass-squared difference some two 
orders of mag nitude greater ~ 1(T 3 - 10" 2 eV 2 0. 

Indeed, the other important new experimental input motivating the present analy- 



sis, is the currently emerging data from the K2K long-baseline experiment [ I3| , tending 
to confirm |T3|] the mass-squared difference from the atmospheric neutrino experiments 
Am 2 ~3x 10~ 3 eV 2 , clearly well above the solar mass-squared difference defined by 
the LMA solar fits and, in particular, subject to the CHOOZ (T§ and PALO- VERDE 
PJ reactor limits on z/ e -mixing, which imply \U e s\ 2 < 0.03, for Am 2 > 10 -3 eV 2 . (U 
denotes the lepton mixing or MNS matrix, fT7||). Note that, in this last respect, the 
new K2K results strongly disfavour the original trimaximal model. 

An obvious solution is to sacrifice the economy of the original trimaximal model 
and acknowledge (in line with most other phenomenological analyses) the existence 
of two distinct mass-squared difference scales, Am 2 >> Am 12 , controlling respectively 
the behaviour of atmospheric and of solar neutrinos Within this context, the totality 
of the data clearly point to a particuluar form for the lepton mixing matrix, which 
turns out to be closely related to the trimaximal hypothesis, and which is given below. 



2. The Trend of the Data 



The atmospheric neutrino results are known to point strongly to twoiold maximal 
Vfj, — v T mixing (or at least to effective twofold maximal — v T mixing at the atmo- 
spheric scale). In particular the SUPER-K zenith angle distribution for multi-GeV 
'//-like' events (Fig. 2a) clearly indicates a suppression of upward with respect to 
downward by a factor of ~ 1/2 (from Fig. 2 the up-to-down ratio for multiGeV 
muons is (U/D)^ ~ 0.53 ± 0.05 for zenith angles |cos#| > 0.20). By contrast, the 
corresponding distribution for 'e-like' events (Fig. 2b) appears to be completely unaf- 
fected, (U/D) e ~ 1.09 ± 0.12. In Fig. 2a the solid curve is the full oscillation curve for 
twofold maximal — > u T mixing for Am 2 = 3 x 10~ 3 eV 2 for a representative neu- 
trino energy, E = 3 GeV. The dashed curve incorporates energy averaging and muon 
angular smearing (with respect to the neutrino direction) and clearly fits the data. 

More generally, the (locally averaged) survival probability < P M > for at in- 
termediate L/E scales, (Am' 2 )" 1 < L/E % (Am 2 ) -1 (where L is the neutrino flight 
path length) is given by the magnitude-squared If/^l 2 of the MNS matrix element 
via < P M > = (1 - |^ 3 | 2 ) 2 + IM 4 , whereby \U^\ 2 = [1 ± ^/l - 2(1- < P M >)]/2 
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(clearly < P M >= 1/2 implies \U^\ 2 = 1/2, as expected). The overall fit to the SK 
atmospheric data (ie. not just < P^ > = (U/D)^ from the multi-GeV muons above) 
gives < P M > = 0.50 ± 0.04 whereby 0.36 < \U^\ 2 < 0.64 (68% CL), certainly consis- 
tent with \Ufja\ 2 = 1/2 and twofold maximal mixing. Independent evidence for strong 
v T mixing comes from the observation of a substantive charged- current v T appearance 
signal, statistically separated from the neutral-current sample in SUPER-K [FS|. 

As indicated in Section 1, intermediate baseline reactor experiments, such as 
CHOOZ and PALO- VERDE in fact provide the best limits on v e (actu- 
ally i>e) mixing at the atmospheric scale (in terms of the vacuum mixing matrix, the 
interpretation of the atmospheric experiments themselves can be seriously obscured by 
terrestrial matter effects, which tend to suppress i/ e -mixing and enhance i/ T -mixing []T9 



in the high-energy limit). Reactor experiments, utilising very low energy (anti) neutri- 
nos and with existing baselines much shorter than the matter-oscillation length in the 



Earth, turn out to be almost completely immune to matter effects PC| . As discussed in 
Section 1, the best reactor limits give |f/ e3 | 2 < 0.03 (95% CL) for Am 2 ~3x 10~ 3 eV 2 , 
consistent with U e s = and thus with two-fold maximal — v T mixing. 

We emphasise that the atmospheric and reactor data do not require U e3 = any 
more than they require (L^l 2 = 1/2 (small non-zero values of U^, and/or somewhat 
different values of |t^3| 2 , eg. \U^\ 2 = 2/3 |2T|, are more-or-less equally acceptable 
experimentally). It is only that U e 3 = and (t^l 2 = 1/2 provide a simple and 
adequate description of the current trend of the data, making twofold maximal — v T 
mixing (for now) the accepted 'default option' |22|| (at the atmospheric scale). 

In a similar spirit, we turn again to the solar data displayed in Fig. 1, drawing 
particular attention now to the solid curve representing the '5/9 — 1/3 — 5/9' bathtub 
(in the LMA solution the base of the bathtub essentially measures \U e 2\ 2 directly). In 
Fig. 1 the data are plotted assuming BP2001 fluxes (§. The SNO , HOMESTAKE and 
SUPER-K data (after correction for the neutral-current contribution in SUPER-K) 
then determine the base of the bathtub, with \U e 2\ 2 ~ 1/3 clearly closely preferred. 
In the flux-independent approach the 8 B suppression is found to be ~ 0.33± g'gy 
(based on the measured SK— SNO difference [l|]). In Fig. 1 the two broken curves 
correspond roughly (in the bathtub region) to the ±1<7 errors on the flux-independent 
supression. Finally, note that with \U e ^ = (or small) the v e survival probability 
outside the bathtub, < P e >, is (inversely) correlated with the value at its base in 
the LMA solution. For vacuum mixing < P e >= (1 — |t/ e 2| 2 ) 2 + |^ e 2| 4 ; so that for 
\U e3 \ 2 = 1/3 we have < P e >= 5/9 ~ 0.56. (Taking account of the 8 B contribution the 
expected gallium suppression is actually closer to ~ 0.53, but this is still consistent 
with the data at the ~ 1.2 a level). Thus the gallium data themselves provide an 
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independent cross-check on the consistency of the LMA solution and on \U e 2\ 2 ~ 1/3. 

As above, we emphasise that the data do not require \U e2 \ 2 = 1/3. If the (implicit) 
energy dependence of the solar suppression is real, certainly |t/ e2 | 2 7^ 1/2 J2"2"| , since 
there can be no MSW bathtub in that case [O. But a somewhat more pronounced 
(eg. a '5/8 — 1/4 — 5/8') bathtub (corresponding to \U e 2\ 2 = 1/4), is clearly far from 
excluded. As before, we regard \U e 2\ 2 = 1/3 as a simple and adequate description of 
the data, which could usefully come to be seen as the default option at the solar scale. 

3. 'Tri-Bimaximal' Mixing 

In this section we simply take the above 'default' values U e3 = 0, \U^\ 2 = 1/2 and 
I U e 2\ 2 = 1/3 as given, and use them to evaluate the resulting lepton mixing matrix. 

The lepton mixing matrix is defined with the rows labelled by the charged-lepton 
mass-eigenstates (e, /1, r) and the columns labelled by the neutrino mass-eigenstates 
(y\, V21 ^3)- Focussing on the last column (the v% column), we note that with U e % = 
and |t^i3| 2 = 1/2, we have |t/ T 3| 2 = 1/2 from unitarity, so that the last column is 
just as in the original bimaximal scheme P5| . Moving to the center column (the v 2 
column), again as a consequence of U e3 = 0, orthogonality requires \U t2 \ = \U^\- With 
\U e 2\ 2 = 1/3 (above) we then have |?7 r 2| 2 = {U^l 2 — 1/3 from unitarity, so that the 
center column is just as in the original trimaximal scheme (see eg. Ref. [[OJ). Finally 
the first column (the v\ column) follows from unitarity applied to the rows. 

Indeed, it was pointed out in Ref. [HJ (even before the SNO data first appeared) 
that a mixing scheme with the z/3 'bimaximally' mixed and the V2 'trimaximally' mixed 
(hence tri-bimaximal mixing) could naturally account for the data, being also discussed 
in the conference literature under the name of 'optimised' bimaximal mixing [p4[| : 
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(where the moduli-squared of the elements are given). The name 'optimised' bimax- 
imal mixing reflected the scheme's pedigree as a special case of the Altarelli-Feruglo 
generalised bimaximal form [25| and its close relationship to the original bimaximal 



form ||23[| . We emphasise that the mixing Eq. 1 is entirely determined by unitarity 
constraints once the above three 'corner' elements are fixed to their default values. 

We should also point out that the mixing Eq. 1 has much in common with the 
Fritzsch-Xing democratic ansatz |21] (which might in fact, see Section 4, reasonably 
be termed 'bi-trimaximal' mixing, as opposed to 'tri-bimaximal' mixing). Indeed the 



5 



Fritzsch-Xing ansatz may be viewed as a permuted form of Eq. 1 (with, somewhat 
remarkably, the crucial prediction U e ^ = made well before the emergence of the 
CHOOZ data \Wl). It should be clear however that the phenomenologies of these two 
mixing shemes are quite distinct, eg. the Fritzsch-Xing ansatz predicted \U e 2\ 2 = 1/2 



and hence no energy dependence [TT1] of the solar suppression, which is now disfavoured 
by SNO The original bimaximal scheme p5| is likewise now disfavoured pUf . 

Asymptotic (L/E — > oo) predictions, specific to tri-bimaximal mixing, are the (vac- 
uum) survival probabilites < P M > = < P T > = 7/18. The corresponding <->• v T 
appearance probability is also 7/18. Note that V e z = implies no Pantaleone reso- 
nance ||27| and no CP violation in neutrino oscillations, which might be considered a 
disappointment experimentally. Nonetheless, it is fair to say that current data 'point 
to Eq. 1, and it is therefore of interest to try to understand what it might imply. In 
the next section we present simple mass-matrices leading to tri-bimaximal mixing. 



4. Simple Mass Matrices 

Our mass-matrices will be taken to be hermitian (ie. we will throughout be implicitly 
referring to hermitian-squares of mass-matrices linking left-handed fields, MM^ = 
M 2 ). Fermion mass-matrices are most naturally considered in a 'weak' basis (i.e. a 
basis which leaves the charged-current weak-interaction diagonal and universal). 

Maximal mixing (whether trimaximal or bimaximal) undeniably suggests permuta- 
tion symmetries P8| . We postulate that, in a particular weak basis, the mass-matrices 
take the following (permuation symmetric) forms: 
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where the real constants a, x, y, z and the complex constants b and b* encode the 
charged-lepton and neutrino masses as follows: 
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with u = exp(i27r/3) and u = exp(—i2ir/3) denoting the complex cube roots of unity. 

In Eq. § the charged-lepton mass-matrix Mf takes the familiar 3x3 circulant form 
[pq| , invariant under cyclic permutations of the three generation indices. The neutrino 
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mass-matrix Mf is real (i.e. symmetric, since our mass-matrices are hermitian) and is 
a 2 x 2 circulant in the 1 — 3 index subset, invariant under the permutation of only two 
out of the three generation indices (generations 1 <-> 3). The neutrino mass-matrix has 
four 'texture zeroes' enforcing the effective block-diagonal form. Note that both 
mass matrices (Eq. |2[) are invariant under the interchange of generation indices 1^3 
performed simultaneously with a complex conjugation. Indeed, in this basis, it is the 
invariance of all the leptonic terms under this combined involution, which guarantees 
no CP- violation here (since Im det [Mf,Mf] changes sign). 

The mass-matrices Mf and Mf are diagonalised by a threefold maximal unitary 
matrix Ui and a twofold maximal unitary matrix U u , respectively: 
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i.e. U}MfUi = diag (m 2 , m 2 



„, m 2 T ) and UlMfU u 



diag (ml, m 2 ,, m 2 ), so that the 



lepton mixing matrix (or MNS matrix) U = UjU u is given by: 
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where the RHS is the tri-bimaximal form (Eq. p]) in a particular phase convention. 
For Dirac neutrinos, the factor of i is readily removed by a simple rephasing of the v% 
mass-eigenstate, yielding tri-bimaximal mixing expressed as an orthogonal matrix: 
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While this concludes our derivation of tri-bimaximal mixing starting from Eq. 0, 
we take this opportunity to remark that one might easily (with perhaps equal a priori 
justification) have interchanged the forms of Mf and Mf in Eq. ^|, taking the neutrino 
mass-matrix to be the 3x3 circulant, and the charged-lepton mass- matrix to be of 
the 2x2 block-diagonal circulant form. Note however that this leads to physically 
distinct mixing which, if the 2x2 circulant is chosen to be in the 1 — 2 index subset, 
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is identically the Fritzsch-Xing democratic ansatz |2T| (hence 'bi-trimaximal' mixing 
as a synonym for the Fritzsch-Xing ansatz, see Section 3 above). We emphasise once 
again that the phenomenologies of the above two mixing schemes are physically dis- 
tinct, with the Fritzsch-Xing ansatz now essentially ruled out, along with many other 
schemes involving energy independent solar solutions |H[ (including trimaximal mix- 
ing 0) following the SNO results 0. 

5. Perspective 

Tri-bimaximal mixing is a specific mixing matrix (Eq. |I]/Eq. |5[) which encapsulates 
the trends of a broad range of experimental data (the LSND oscillation signal 



was not considered on the grounds that it still awaits confirmation). Tri-bimaximal 
mixing is closely related to a number of previously suggested lepton mixing schemes, 
notably trimaximal mixing ||, bimaximal mixing |23[], the Fritzsch-Xing democratic 
ansatz |^TJ and the Altarelli-Feruglio scheme |25| (of which tri-bimaximal mixing may 
be considered a special case). For the future, the couplings of the heavy neutrino, z/ 3 , 
are expected to be measured more precisely in long-baseline experiments like MINOS 
[ P^l (and other projects Pl|). In particular, the limits on \U e3 \ 2 should continue to 



improve. Regarding the v 2 couplings, if tri-bimaximal mixing is right, the KAMLAND 
experiment ]35| should confirm the LMA solution, measuring a v e survival probability 
tending to < P e >= 5/9 ~ 0.56 at sufficiently low-energy. Corresponding predictions 
for Up disappearance and — > v T appearance, see Section 3, look very hard to test 
(experiments with necessitate higher energies and hence longer baselines, with mat- 
ter effects generally dominant over vacuum effects for V2). Finally, exact tri-bimaximal 
mixing would imply no high-energy matter resonance and no (intrinsic) CP-violation 
in neutrino oscillations, which might be considered a disappointment experimentally. 
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Figure Captions 



Figure 1. The solar data, including the recent SNO measurements, plotted versus 
neutrino energy E. The gallium points are plotted at < 1/E > _1 ~ 0.5 MeV, the 
SUPER-K point at < 1/E > _1 ~ 10 MeV and the HOMESTAKE point (which in- 
cludes an ~ 15% contribution from the 7 Be-line at ~ 0.9 MeV) is plotted (somewhat 
arbitrarily) at < 1/E > _1 = 5 MeV. The solid curve is the expected energy depen- 
dence for tri-bimaximal mixing, Eq. 1, with Am' 2 = 5 x 10 -5 eV 2 , which predicts a 
suppression of 1/3 in the base of the bathtub and 5/9 outside. 

Figure 2. Atmospheric neutrino data. The multi-GeV zenith angle distributions for 
a) /i-like and b) e-like events in SUPER-K. The solid curve is the full oscillation curve 
for tri-bimaximal mixing (Eq. 1/Eq. 6) with Am 2 = 3 x 10 -3 eV 2 for a representa- 
tive neutrino energy E = 3 GeV, and the dashed curve shows the effect of angular 
smearing with respect to the neutrino direction and averaging over neutrino energies. 
Tri-bimaximal mixing (or indeed any mixing hypothesis which is effectively twofold 
maximal — v T mixing at the atmospheric scale) fits the data well. 
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